Particle dynamics in a relativistic invariant stochastic medium. 
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The dynamics of particles moving in a medium denned by its relativistically invariant stochastic 
properties is investigated. For this aim, the force exerted on the particles by the medium is defined 
by a stationary random variable as a function of the proper time of the particles. The equations 
of motion for a single one-dimensional particle are obtained and numerically solved. A conserva- 
tion law for the drift momentum of the particle during its random motion is shown. Moreover, 
the conservation of the mean value of the total linear momentum for two particles repelling each 
other according with the Coulomb interaction is also following. Therefore, the results indicate the 
realization of a kind of stochastic Noether theorem in the system under study. Possible applications 
to the stochastic representation of Quantum Mechanics are advanced. 
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I. INTRODUCTION 

Stochastic processes made their appearance in research 
in Physics long time ago and their theory has played an 
important role in the description of systems which do not 
behave in a deterministic manner 12 3 . In particular, the 
study of the dynamics of particles lying inside material 
media has been the object of high interest. A classi- 
cal example is the study of the Brownian motioni. A 
large deal of those investigations had a non-relativistic 
character and the random interactions with the back- 
ground medium were considered as being dependent of 
the state of motion of the particle, that is, lacking in- 
variance under the changes of the reference systemic. 
Another large class of studies in this field had been di- 
rected to show the equivalence with random processes of 
the solutions of quantum relativistic or non-relativistic 
equations, like the Klein-Gordon, Dirac and Schrodinger 
onesi 4 i 5 i 6 i 7 i 8 i 9 i 10 ' 11 i 12 . Two basic additional subjects in 
connection with stochastic processes in Quantum The- 
ory are: The attempts to derive the 'collapse' of the wave 
function during measurements from the existence of ran- 
dom perturbations in Quantum Mechanics (QM)^ 4 *^^* 1 ^, 
and the study of the decoherence processes and their role 
in spontaneous transitions from pure to mixed statesit 

The main objective of the present work is to inves- 
tigate some consequences on the motion of a particle 
determined by the action exerted over it by a medium 
which random properties are defined in a relativistically 
invariant form. The basic motivation is simple: It is 
recognized that the Copenhagen interpretation of Quan- 
tum Mechanics (QM), is the most successful and dom- 
inant from all the existing ones. However, it is also ac- 
cepted that its approach to measurements constitute one 
its more puzzling aspects, which up to now is widely de- 
bated in the literaturo 4 ! 5 ! 6 ! 7 ! 8 ! 9 ! 10 ! 11 ! 12 ! 13 ! 14 ! 1 -^ 16 . Let us 
suppose for a moment, that in opposition to the Copen- 
hagen interpretation and in accordance with Einstein ex- 
pectations, the particles in Nature are in fact localized 



at definite points of the space at any moment. Then, 
the only way we can imagine for the quantum mechan- 
ical properties of the motion to emerge from a model, 
is that the action of the vacuum on the particles have a 
stochastic character. But, the relativistic invariance of 
the vacuum, leads to expect that the acceleration felt by 
the particle in its proper frame should be a stationary 
random variable as a function of the proper time. This 
circumstance motivates the study started here about the 
motion of particles inside a random media showing the 
above mentioned property. For the sake of simplicity 
the one dimensional motion is considered. It is not at- 
tempted to show the equivalence of the dynamics in the 
medium with the one predicted by the quantum mechani- 
cal equations. The purpose in this first step, being redun- 
dant, is to study general properties of the motion of one 
and two particles assuming two main conditions: a) They 
have a definite localization in the space at any moment 
and b) The forces acting on the particles have the above 
mentioned random properties which are independent the 
observer's inertial reference frame. 

The work will proceed as follows. Firstly, the equations 
of motion of the particles under the action of the medium 
are formulated. For this purpose the properties which 
ensure the relativistic invariance of the motion under the 
action of the medium are stated by specifying the form of 
the random forces. Further, the equations of motion of 
a single particle are written and solved and a statistical 
analysis of the random properties is done. A main 
conclusion following is the existence of a conservation 
law for a mean drift momentum and kinetic energy of a 
'free' particle propagating in the medium. It indicates 
the validity of a kind of stochastic Noether theorem which 
links the relativist invariance of the stochastic motion 
with the conservation of the mean 4-momentum of the 
particle. 

Further, the conservation law is studied for the 
mean of the addition of two four momenta associated to 
the scattering of two identical particles, which repel each 
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other through an instantaneous Coulomb interaction. 
It is concluded that the random action of the medium 
does not modify the usual conservation law, valid for the 
impact in the absence of external forces. 

A review of the results and future extensions of the 
work are presented in a conclusion section. Some pos- 
sibilities to extend the study are advanced. In general 
terms, our view is that the form of the analysis have the 
chance of being useful in the search for consistent hidden 
variables models. The study of these possibilities will be 
considered elsewhere. 



II. EQUATION OF MOTION 

In this section we will obtain and solve the Newton 
equation of motion for a particle on which a random force 
F p (r) is acting. A one dimensional system will be con- 
sidered to make the discussion as simple as possible. The 
force will be defined as a vector in the proper reference 
frame of the particle and will depend on the proper time 
t. That means, in each instant we will consider an in- 
ertial reference frame moving relative to the observer's 
fixed frame (Lab frame) with the velocity of the particle 
v and having the origin of coordinates coinciding with it. 
In this system of reference, after a time increment dr, it 
is possible to write 

F p (t) dr = mo dv', (1) 

where too is the proper mass of the particle. 

The particle reaches a small velocity dv' relative to 
this system and a new velocity respect to the Lab frame 
v + dv, given by the equation 

-dv 1 



dv = 



1 + 



{v + dv') X (1 



However, since the relativistic invariance condition will 
be imposed on F p (r) this will be the type offeree mostly 
considered in what follows. Integrating the equation[3in 
the proper time it follows 



F p (t) dr + C — toq 



1 



dv, 



which determines the velocity the Lab frame v(t) as a 
function of the proper time r, only through the depen- 
dence of r, of the integral of the random force F p (t). The 
explicit form of v becomes 

v(t) = c ■ tanh[— • ( / FJt) ■ dr + &)], (8) 

TO C J 



v dv' 
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where C is an arbitrary constant. 



III. THE RANDOM FORCE 

As it was mentioned in the Introduction, the medium 
under study will be defined in the proper frame as ran- 
domly acting over the particle being at rest in it. That 
is, its action in this reference system will be given by a 
stochastic process showing no preferential spatial direc- 
tion and assumed to be produced by an external relativis- 
tic system which dynamics is unaffected by the presence 
of the particle. Its is also natural to impose the coin- 
cidence of the distribution function of the forces of the 
medium for a large sampling interval of proper time T 
and the one obtained fixing the proper time r, produced 
by an ensemble of a large number of samples of the forces 
taken during long time intervals T. 
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where c is the velocity of light, 
speed in the Lab frame dv is 

2 



(3) 

Thus, the variation of 



dv'. 



(4) 



From expressions (QJ and (0J the required differential 
equation for the motion is obtained : 

dv 



F p (t) 



to 



d x(t) 



(1 



d x 
] J-r 



(5) 



(6) 



d t 

It is useful to state the relation between the strength 
of the force in the Lab system and its proper frame coun- 
terpart, which is: 



F P (r) 



1 - -F l (t). 

c 2 - 



(7) 




FIG. 1: The figure shows a realization of the force field 
corresponding to a spectrum of N = 200 frequencies. The 
horizontal coordinate is the time in seconds and the verti- 
cal one is the force in Newtons. The amplitude was fixed 
fo = 1 (N t ) (N t = Newton). 
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These conditions, can be assured by a random force 
F p (t) being stationary, ergodic and symmetrical dis- 
tributed about the zero value of the force. A numeri- 
cal realization of a band limited white noise distribution 
obeying these properties is implemented in reference^ 
and will be employed here. Concretely, the expression for 
the stochastic force given in the proper reference frame 
as a function of the proper time will be taken in the form 

fa N 

/jv(t, (fi) = — ^cosOi (N) T + <pi), (9) 

Wi(N) = ^, i=l,...N, (10) 

where the N phases (fi, i = 1, ...N, are randomly chosen 
with a uniform distribution between and 2tt. The in- 
teger number N is the number of frequency components 
of the numerical band limited white noise distribution. 
The bandwidth will be chosen as a fixed one and equal 
to Aw = 8ir. Clearly, the exact randomness for an 
arbitrarily large time interval only will attained in the 
infinite limit of N. The parameter /o controls the ampli- 
tude of the force values. Note that the absence of a zero 
frequency component is implied by the condition of the 
process not showing a preferential direction in space. 
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This result is natural due to the fact that the force is 
normalized and its absolute value can not exceed /o. 

IV. A PARTICLE IN THE MEDIUM 

In this Section the existence of a conservation law for 
the mean momentum (to be also called from now on the 
drift momentum) will be shown for a particle moving in 
the before defined random medium. To this purpose 
let us employ the solution © which was found in the 
previous section We will combine this expression with 
the results obtained from the definition of the random 
force ensemble in equation @, to explicitly determine 
v = v(t). The resulting relation links the velocity of the 
particle in the Lab frame with the proper time measured 
by a clock fixed to it. Once having v(t) we will comment 
about its stationary random behavior. The existence of a 
non vanishing conserved mean value of the velocity < v > 
will be shown. Starting form the method defined above, 
the velocity distribution function will be determined for 
the frequency spectrum defined before and few represen- 
tative values of the arbitrary constant C . After that, the 
position of the particle x(t) as a function of the time t 
measured at the Lab frame will be evaluated. 




FIG. 2: The distribution function of the ensemble of forces 
for the set of frequencies w, = Jgjji (s _1 ), i = 1, ...250 and 
fa = 1 (Nt). 

Figure ^ shows the force distribution fx for a 
realization of the ensemble of forces for N = 200. 
A realization here is called the time evolution of the 
force, for the set of randomly fixed phases ipi at the start. 
The picture qualitatively shows the stationary charac- 
ter of the random force. Figure |3 depicts an interpola- 
tion curve of the data for the distribution function w(f) 
corresponding to an ensemble of random realizations 
of the force. The white noise frequency spectrum was 
defined by N — 250 frequency components Wi — 
(s- 1 ), % = 1,...250. 

The force amplitude fixed was as fo = 1 (Nt). The 
picture corresponds to a large sampling time T ( it is 
sufficient to be at most equal to the period of the smallest 
frequency of the spectrum) . Notice the even character of 
the distribution and its rapid decay to zero. Of course, 
this occurs inside the interval defined by /1.2 = ±1 Af(.. 



FIG. 3: The velocity of the particle in the laboratory frame 
against the proper time and the arbitrary constant C. The 
parameters for the random force were S J^ - = 0.1 (s" 1 ), 

Wi = S TT (s- 1 ), i = l,..-,250. 

Taking F p (t) to be given by the white noise force 
/n(t) we have : 

N 

F P (r) - Mr) = ^ Yl C0S ^ T + V*)" ( U ) 
Then, integrating with respect to r produces 

If p (t)= f F p (t) dr = ^ f [£cos(^T + ^)]dT, 
Jo Jo i=1 
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Substituting l(T2*jl into (JSJ, gives 



u(t)=c tanh[ (J^(r) + C)], 

moc 

iV 



(13) 



= c tanh[— (A[^Isin(^r + ^)]+C)]. 



(14) 
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FIG. 4: The velocities of the particles (divided by c) in the 
laboratory system vs. the proper time for four specific values 
of C: a) (7=0 b) (5=0.3 c) C=-0.2 and d) C=-0.8. Note 
the non vanishing value of the mean velocity for C different 
from zero. 



It can be seen that the summation within the argument 
of the hyperbolic tangent is symmetrically distributed 
around its zero value. Since the tanh(x) is an anti- 
symmetric function, it follows that when C / there 
will be a nonvanishing mean velocity of the particle in 
the medium. As the mean value of the kinetic energy is 
also conserved, it can be said that the mean 4-momentum 
of the particle is conserved. Moreover, the relativistic 
invariance of the system also implies that the mean 4- 
momentum of the solutions for any two values of C should 
be linked by certain Lorentz transformation. Thus, the 
particle trajectories for the various values of C are sim- 
ply a fixed trajectory (in the stochastic sense) after being 
Lorentz transformed into certain moving frame. 

The picture in figure shows the behavior of 
v = v(t, C) for the following values of the parameters 



Jo 



0.1s- 



w. 



8tt i 



1,...,250. Note 



Stoi c > — i N ' 

that for fixed values of C (that means, for certain ini- 
tial conditions) the velocity rapidly oscillates around a 
value being close to the quantity tanh[(7]. It illustrates 
the mentioned conservation of the mean drift velocity 
for the particle in spite of the random oscillations of the 
medium. Figure 01 shows the same dependence for spe- 
cific values of C . 
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FIG. 5: The distribution functions: a) w(v) for (7=0.1, b) 
w(v) for (7=0.3 and c) w(v), for (7=0.5. Note the distortion 
of the symmetry around the center, when C grows. 



A. Distributions 

Let us consider now the numerical evaluation of the 
distribution function w(v) (for measuring a given value 
of the velocity v) for few representative values of C . This 
function is found on the basis of the ergodic property of 
the system, that is, by performing a large number of 
evaluations of the velocity's expression with time run- 
ning from zero to a "sampling" value T, for afterwards 
compute the number of times, for which v takes val- 
ues in a small neighborhood of a given value. Further, 
after interpolating the results, the distribution functions 
are obtained after dividing by the fixed size of the men- 
tioned neighborhood. In figure (JjjJ the distribution w(v) 
is plotted for a few values of C. 

These pictures correspond to the frequency spectrum 
Wi = i = 1--.250, but with the amplitude fixed 

by 4 * 25 = I- The almost independence of the form of 
these curves on the size of a large "sampling" time T 
employed (whenever N is sufficiently large) indicates the 
approximate validity of the ergodic property of the white 
noise implementation employed. 

Note the rapid decay away from the mean value and 
the deformation of the symmetry around it, of the dis- 
tribution w(v) for different values of C. These distribu- 
tions allow to calculate the mean velocities in each of 
the cases. This method will be employed in the next 
sections for finding the mean values as integrals over the 
domain of the quantity, of its value times the distribution 
function. 



B. Particle velocities and trajectories 

It can be noted that the explicit solution for the veloc- 
ity obtained in JSJ) corresponds to this velocity in terms 
of the proper time r. However, in order to integrate 
the velocity to find the random particle trajectories as 
functions of the time in the laboratory frame, it becomes 
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necessary to know the functional relationship r(t, C) for 
the considered trajectory. Then, let us consider now the 
numerical determination of this relation. The parameters 
of the random force will be s J° 1qC = 0.1 (s _1 ) , Wi = ^ i 

(s _1 ), i = 1, 2, ...,250. Using the solution v = v{t, C) it 
follows 

V c ' 

t(0) = 0, (16) 

and employing it, the values of t(r, C) were numeri- 
cally found. Afterwards, finding the inverse mapping 
t = r(t,C), this function is depicted in figure with 
C as a parameter. Finding the composed function for 



the values of the positions x(t) respect to the laboratory 
reference frame are obtained. The trajectories of the 
particles are shown in figure |H1 for the same set of values 
of C. 
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FIG. 6: Dependence on t of the proper time in the labora- 
tory r = r(t,C), with (7 running along a third axis. Note 
that, since the dependence, by definition, always should be 
monotonous, the random oscillation are not apparent. They 
are however present in the local slope of the curves. 

various values of C , the observer's time dependence of 
the particle velocity u(t, C) follows. The results are 
illustrated in figure [7| for the chosen values of C. 
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FIG. 7: The velocities in the laboratory frame (divided by c), 
but now plotted as functions of the time t in this same frame 
of reference, for the values of (7: a) (7 = 0.01, b) (7=0.4, c) 
(7=-0.03, d) (7=-0.5. 

Further, after integrating the calculated velocities with 
respect to t (starting at t = and imposing x(0) = 0) 



FIG. 8: The particle coordinates divided by c: - = -x(t) as 
functions of the time in the laboratory frame. In a) and b) C 
= 0.01 but the time scale is different. For c) and d) (7=0.4 
and again the two scales are different. Note that for C small 
the randomness seems to be greater than for larger (7. 

Note the randomness of the motion in the case of small 
absolute values of C which is not the case for the larger 
ones. This property can be understood analyzing the 
pictures in figure JJJ). In the cases in which the veloc- 
ity constantly changes from positive to negative values 
and vice versa, the randomness is more evident. Thus, 
the aleatory effect is more visible in one case than in the 
other, because the curve changes from a monotonous to 
non monotonous behavior. But, after taking into ac- 
count the relativistic invariance of the model, as noticed 
before, it can be seen that the curves for different values 
of C should transform any one into another by a Lorentz 
transformation (exactly, when the same realization of the 
force is employed, or statistically, if another realization 
is used). Therefore, the apparent increased randomness 
for low C is only a visual effect enforced by the change 
of sign of the velocity in frames in which its mean value 
is sufficiently low. 

V. TWO PARTICLES IN THE MEDIUM 

Let us consider in this final section the conservation 
of the mean total momentum of a 'closed' system of two 
particles which travel in the medium by also interacting 
between them. A repulsive interaction between the two 
identical particles 1 and 2 will introduced. The forces 
between the particles will be defined in the laboratory 
frame for studying the impact between them there. They 
will be chosen as satisfying the Law of Action and Reac- 
tion and having the Coulomb form. If, for example, the 
forces have an electromagnetic origin, the retardation ef- 
fects will be disregarded in order to assure that the field 
will not deliver momentum to the set of two particles. 
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This approximation is appropriate for low particle veloc- 
ities in its random motion in the Lab frame. 

The system of differential equations will be written and 
numerically solved. The results for the position and ve- 
locity of both particles x\ (t) , v\ (t) , x 2 (t) , v 2 (t) , will show 
how after the impact, any of them deliver to the other 
its mean drift velocity and its type of randomness. As 
it was mentioned above, the expression for the repulsive 
force in the Lab frame will have the Coulomb form 

F repi {x ll x 2 ) = ~F rep2 (x 1 ,x 2 ) = | — ■ — 13(2?!- 2? 2 )- 

\Xi X 2 \ 

(17) 



Consider now the relativistic Newton equations for 
both particles in their respective proper frames and also 
the two transformations between the common laboratory 
time t and the two different proper times t\ and r 2 . 
These relations may be written as 



I 1W = F ^ + ^ - X -^-) i F rePl (xi ) x 2 ) ) (18) 

r^I ^ = F pM + (!- -^)* W^), (19) 
- (1 ~ n(t ) = t 1o , xi(t ) = x lo , x'xito) = vi , (20) 



dt 

= (1 _ 



(l-^^) 2 , 7a(*d) = 7a 0) aj 2 («o) = isao, a/a(*o) = «a - ( 21 ) 



In order to be consistent with the non retarded approx- 
imation for the Coulomb repulsion, as mentioned above, 
initial velocities for the particles being relatively small 
with respect to c were taken. This assumption makes 
the considered approximation to be satisfied if the inten- 
sity of the repulsive force is weak so that the velocities 
after the impact will also be small. Employing the 
white noise parameters Wi — §§ji, i = 1,...250 and 
f° 125 = 1, we obtained the numerical solutions shown in 
figure El for various initial conditions. The conservation 
law of the total mean momentum of the two particles 
can be noticed. Figure ED shows the time variation of 
the velocities of the particles. The process of exchange 
of their drift momenta is clearly illustrated and therefore 
the corresponding conservation of the total momentum. 
This outcome follows in spite of the strong oscillations 
produced by the action of the medium, in opposition to 
what happens within a standard material media. In these 
systems the stochastic action normally tends to stop the 
motion of the particle making the drift velocity to vanish, 
assumed the absence of accelerating external fields. The 
calculated values for the mean velocities before and after 
the shock are 

< v Xo > = 0, < v 2o >= 0.152, (22) 
<v lf > = 0.152, < v 2f >= 0, (23) 

which numerically confirms the conservation of the to- 
tal momentum in the shock of two particles forming a 



'closed' system immersed and moving in a relativistic 
random medium. 



VI. CONCLUSIONS 

The one-dimensional equation of motion for a particle 
moving in a medium having relativistic invariant stochas- 
tic properties is formulated and numerically solved. The 
velocity of the particles in the medium is a function of the 
proper time only through the integral of the force in the 
proper reference frame. This relation directly shows the 
existence of a stochastic conservation law: A free par- 
ticle in the defined random medium conserves its mean 
momentum and kinetic energy along its motion. 

The problem of two particles moving in the medium 
is also investigated after considering a momentum con- 
serving Coulomb repulsion between them. The evalu- 
ated solutions generalize the conservation law: The sum 
of the two drift momenta of two particles moving in the 
medium without any other external action conserves after 
a shock. 

An interesting outcome is that for two different shocks 
with identical initial conditions for both particles when 
they are far apart, the drift velocities before the impact 
are not uniquely determined by the initial conditions. 
These velocities also depend on the array of phases uti- 
lized for the realization of the force. Therefore, it follows 
that the results of the impact will show a dependence 
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FIG. 9: The trajectories of the particles xi(t) and X2(t) 
during the impact for the repulsion parameter a = 0.01 (N x 
m 2 ) and initial conditions m(500) = £2(500) = 0, «i(500) = 
0.015, t> 2 (500) = 0.15, n(500) = 0, r 2 (500) = 0. The curve 
at the bottom is the same that the one at the top, but at a 
lower scale in the vertical coordinate. 
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FIG. 10: The time dependence for the velocities (divided by 
c) of both particles: a) vi (t) jc and b) V2 (t) /cm the considered 
shock. Note, the exchange of their mean velocity values. 

on the specific realization of the random force. This cir- 
cumstance implies that the probability distribution as- 
sociated to an ensemble of particles all situated at the 
beginning at the same point and with a fixed value of 
the velocity, will not define only one direction of the 
trajectory for large times. The probability distributions 
of such an ensemble, on the other hand, should evolve 
in space and time in a rclativistically invariant manner, 
as the Lorentz invariance of the system indicates. The 
same conclusion can be traced for any other sort of fixed 
boundary conditions. Therefore, the set of possible ini- 



tial conditions for the particle (considering also that they 
can be placed in the medium at different initial spacial 
and temporal points for the construction of the ensemble) 
should be expected to be equivalent to the set of all pos- 
sible space time evolving ensembles that can be observed. 
The above remarks suggest some possible extensions of 
the present work, which are described below: 

• The indicated dependence of the ending results of 
the shocks, not only of the initial conditions, but 
also of the concrete realization of the random force, 
led us to think in extending the results to the 2D 
and 3D cases. The idea is to study the spatial and 
temporal behavior of the ensemble of outcomes o a 
series of shocks "prepared " with fixed initial mean 
velocities for both particles, and to compare the 
results with the corresponding predictions of the 
quantum scattering. 

• Another task of interest which is suggested by this 
study is to investigate the existence of preferential 
bounded states in the case that the Coulomb in- 
teraction is considered as attractive. A particular 
simpler situation could be to assume one of the par- 
ticles as very massive, that is, merely acting as an 
attracting center. In both cases the study could 
consist in determining the statistical properties of 
the stationary trajectories given the initial condi- 
tion. 



In conclusion, we would like to underline that it seems 
not without sense that the realization of the above pro- 
posed studies could be of use in the justification or search 
for hidden variable approaches to Quantum Mechanics 
(QM). As it can be seen from the discussion, the resulting 
models could not show the limitations of the Bohm ap- 
proach (like the absence of predictions for the " guided" 
motions of the particles for real wave functions, for exam- 
ple). Moreover, the analysis seems of interest in seeking 
for a theory not requiring to fix a particular random pro- 
cess to each stationary state, but one in which all the sta- 
tistical properties can emerge from a general framework. 
In the imagined outcome, the particle could propagate 
having a probability for transit or stay into each one of 
a set of stationary random motions which could be asso- 
ciated (but now within the general context) to the par- 
ticular eigenfunctions of the system. As for the allowance 
of the necessary generalizations needed to make contact 
with reality, it can guessed that the possible generaliza- 
tion of the statistical Noether properties (detected in the 
simple model considered here) could lead to the conserva- 
tion of the mean values for the angular momentum and 
other internal quantities, in analogy with the case of the 
linear momentum. We think these possibilities deserve 
further examination that will be published elsewhere. 
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